ON BUDAEV AND BOGY'S APPROACH TO DIFFRACTION BY THE 2D 
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Abstract . Several semianalytical approaches are now available for describing diffraction of a plane wave by the 2D 
(two dimensional) traction free isotropic elastic wedge. In this paper we follow Budaev and Bogy who reformulated the 
original diffraction problem as a singular integral one. This comprises two algebraic and two singular integral equations. 
Each integral equation involves two unknowns, a function and a constant. We discuss the underlying integral operators 
and develop a new semianalytical scheme for solving the integral equations. We investigate the properties of the 
obtained solution and argue that it is the solution of the original diffraction problem. We describe a comprehensive 
code verification and validation programme. 
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1. Introduction. Sommerfeld (1896) was the first to solve a wedge diffraction problem - that of 
diffraction of an electro-magnetic wave by the perfectly conducting semi-infinite screen. In this famous 
paper, he obtained an exact solution in the form of a superposition of plane waves propagating in 
complex directions. Nowadays this representation is called the Sommerfeld integral. Sommerfeld 
(1901) also solved the problem of diffraction of an electro-magnetic wave by the wedge whose angle 
is a rational multiple of 180°. He suggested that since any irrational number can be represented as 
the limit of a rational sequence the result could be extended to any wedge angle. The procedure was 
implemented by Carslaw (1920). 

Another analytical method for solving two-dimensional diffraction problems was developed by 
Smirnoff and Sobolev (1932). They tackled the impulse diffraction by an acoustic wedge under Dirich- 
^r^ . let or Neumann boundary conditions, without recourse to the frequency domain. Chapter 12 in Frank 

and von Mises (1937) which has been written by Sobolev, contains a detailed exposition of the method 
and describes its application to the acoustic wedge. Friedman (1949a, 19496) and Filippov (1959) both 
applied this approach to modelling diffraction of an elastic wave by the linear semi-infinite crack in 
an elastic plane. The time-harmonic version of the problem was treated by Maue (1953) with the 
r^ ' Wiener-Hopf technique. 

In the 1950's Malyuzhinets (1955-1959) solved the problem of acoustic plane wave diffraction by 
the wedge with impedance boundary conditions. He sought a solution in the form of the Sommerfeld 
integral - whatever the wedge angle. Then the boundary conditions have been reformulated in the 
form of functional, difference, equations F(u ± a) — R (uj)F(— u> ± a), where F is an unknown 
meromorphic function and R = (—sin u — a ± )/(sin u> — b ), with a and b known constants. 
Malyuzhinets has shown that these equations have an analytical solution. Williams (1959) solved the 
same problem independently: he took advantage of the structure of R 1 * 1 to simplify the Wiener-Hopf 
factorisation. 

A good exposition of Malyuzhinets' theory as applied to the acoustic and electro-magnetic wedges 
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was made by Osipov and Norris (1999). The authors have emphasised the features of the theory which 
pertain to impedance, the higher order, boundary conditions. These lead to a more involved form of 
R which nevertheless remain rational functions of sin u> and cos u>, allowing for an analytical solution 
(see also Tuzhilin 1973.) 

Several mathematicians worked on elastic wedge problems throughout the 1960's and 70's. An 
analytical solution was found for the slippery rigid elastic wedge - the wedge with the zero tangential 
traction and zero normal displacement on the boundary (e.g. Kostrov 1966). The three-dimensional 
smooth elastic wedge with mixed boundary conditions also was treated analytically (e.g. Poruchikov 
1986). A comprehensive review of early papers dealing with diffraction by the elastic wedge was made 
by Knopoff (1969). 

For the traction-free elastic wedge an analytical solution has been found only in the degenerate 
cases of wedge angles of 360° (a linear semi-infinite crack - see Friedman 1949a, 19496; Filippov 1959 
and Maue 1953) and 180° (an elastic half- plane). Plane shear wave incidence along the bisectrix of a 
quarter plane can also be treated analytically. 

There have been many attempts to find an analytical description of diffraction by the elastic wedge 
in non-degenerate cases. The problem became a diffractionist's analogue of the famous Fermat's Last 
Theorem. " Diffraction-fermatists" have not been able to obtain any promising results. It appears 
that there can be no analytical solution. The traction-free elastic wedge problem has to be tackled 
numerically. 

Since the mid 1980's many researchers followed this route. Several relied on potential theory to 
reduce the problem to integral equations. In particular, the problem of Rayleigh wave diffraction by 
the traction-free elastic wedge was attacked this way by Gautesen (1985 - 2002c) and Fujii (1994 and 
references therein). Fujii (1994) studied the wedge angles between 36° and 180°. Gautesen started with 
the wedge angles of 90° (Gautesen 1985, 1986 and 2002a), 270° (Gautesen 20026) and then considered 
the wedge angles restricted to the interval 189° to 327° (Gautesen 2001) and to the interval 63° to 
180° (Gautesen 2002c). A detailed exposition of a version of the boundary integral equation approach 
called the spectral functions method has been given by Croisille and Lebeau (1999). The authors 
consider the challenging problem of diffraction of a plane acoustic wave by the elastic wedge immersed 
in liquid. The monograph contains both fundamental developments and numerical results. The ideas 
of Croisille and Lebeau have allowed Kamotski and Lebeau (2006) to reformulate radiation conditions 
that are consistent with physical considerations and also allow one to prove theorems of existence and 
uniqueness. 

In another series of recent papers the problem of diffraction by the traction-free elastic wedge was 
reduced to the singular integral equations by representing the elastodynamic potentials in the form of 
the Sommerfeld integrals - in the spirit of Malyuzhinets' approach. Larsen (1981) appears to have been 
the first to attempt to apply this approach to describe diffraction by the elastic wedge. He considered 
the zero displacement boundary condition and reduced the problem to a system of functional equations 
in analytic functions. He suggested that Chebyshev's polynomials could be used to solve the problem 
numerically, but published no further results. Budaev (1995) and Budaev and Bogy (1995, 1996 
and 2002) have taken the Sommerfeld-Malyuzhinets approach further and reduced the problem to 
a system of two functional equations in two analytic functions. These equations are similar to the 
ones obtained in the impedance acoustic wedge problem mentioned above (Malyuzhinets 1955-1959) 
but are matrix rather than scalar. Budaev and Bogy (1995-2002) presented numerical results on 
scatter of an incident Rayleigh wave, but explanation of their numerical scheme and some theoretical 
arguments are vague. Their computed reflection and transmission coefficients do not always agree 
with the corresponding numerical results obtained by other authors. Nevertheless, as we show in this 
paper, Budaev and Bogy's approach is valid. Our immediate aim is to clarify certain aspects of their 
theoretical treatment, reduce the wedge problem to a singular integral one and analyze its properties, 
develop a stable numerical procedure for its solution and then, assuming the incident field Rayleigh, 
evaluate the corresponding diffraction, reflection and transmission coefficients. We argue that the 
obtained solution is that of the original problem. 

In §Sj5|and|21we outline our own semianalytical recipe for solution of the singular integral problem, 
and in 21we verify and validate the resulting code. In Appendix \j\ we describe the nomenclature and 
in other Appendices, offer the necessary theoretical considerations, formulas and numerical options. 
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Figure 1.1. Geometry of the traction free elastic wedge. 
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Figure 2.1. The Sommerfeld contours in the complex u-plane. 



2. Statement of the problem and the Sommerfeld amplitudes. Let us briefly present the 
full statement of the original diffraction problem. We seek the elastodynamic potentials (pi — <fii (kr, 9) 
that satisfy the Hclmholtz equations in the two-dimensional the wedge of angle 2a with traction free 
faces, that is, we address the boundary value problem 
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Above and everywhere below, the parameter k is the shear wave number; 7 = cs/cp is the ratio of the 
shear and compressional speeds cs and cp and the subscript i takes values or 1. The geometry of 
the problem is shown in Fig. 11.11 Given an incident wave we seek the scattered potentials satisfying 
the radiation conditions at infinity (analogous to the ones in Kamotski and Lebeau 2006, Theorem 
4.1) and bounded elastic energy condition at the wedge tip. 

Note that the potentials are related to displacement u = u(x) via u = V0o + V ^1, where the 



nabla operators are V = (8 X 
Hclmholtz equations imply 



8 y ), V = (8 y , —8 X ). Note too that in view of this representation the 



<Mx) = -( 7 fcr 2 V-u(x), 0i(x) = -fr 



-2 V _L 



u(x). 



(2.4) 



It follows that 4>i(kr,9) are uniquely defined by u(x) and since for the corresponding Lame problem 
in u(x), the existence and uniqueness results have been proven (Kamotski and Lebeau 2006 and 
Kamotski 2003, Theorem 3.1), the above diffraction problem as formulated in terms of the potentials 
has a unique solution too. Moreover, due to (|2.4J) all the analytical properties of u(x) established in 
Kamotski and Lebeau (2006) imply the corresponding properties of c/>i(kr,9). 

The solutions (f>i (kr, 9) of the Hclmholtz equations can be represented in the form of the Sommer- 
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feld integrals 



<j>i{kr,0)= / <P l (cj + 8)e i " /krcoB "dcj, (2.5) 

cue 



which can be rewritten as 



4>i(kr,6)= $i(u + 0) - $i(-uj + 6) c^ krcos "diu. (2.6) 

J c L J 

Here the II shaped contour C runs from — ir/2 + ioo to 3tt/2 + ioo. The contour C is the reflection of C 
with respect to the origin, and the full integration contour C(J C is shown in Fig. 12.11 One can justify 
(|2.6(l by applying the Laplace transform in kr to <fii(kr, 9), changing the Laplace variable stow, such 
that s = i cos u and exploiting the behavior of 4>i{kr,&) both in the vicinity of the wedge tip and 
at infinity (for details see Appendix EJ. It follows that the analytic functions &i(uj) possess a finite 
number of singularities in any vertical strip and are regular at the imaginary infinity. Moreover, in 
Appendix [B] we show that the tip behavior of 4>i(kr, 9) determines the asymptotic expansion of $»(w) 
at the imaginary infinity and provides us with an asymptotic estimate 

*i(w) =0(e- Ro p Im "I), Rep>-1, |Imw|->oo, (2.7) 

which assures the convergence of integrals (|2.6H . Note that these integrals are invariant under the 
transformation $j(w) — > &i(w) + const. In Appendix iBl we also show that as Im u> — > oo, the 
asymptotic expansion of the preexponential factor in (|2.6|l contains a constant term. Below, we choose 
so that the constant terms in their asymptotic expansions valid as Im lo — > oo and Im lo — > — oo differ 
only by sign. We call those $i(w) the Sommerfeld amplitudes. 

Introducing the closed contour C(J(— C2) U C[j(— Ci), the Sommerfeld integrals l|2.5[) can be 
evaluated as the contributions of the poles and branch points inside that contour plus the sum of 
integrals over the steepest descent contours Ci and C2. The former describe all bulk and surface 
waves arising in the problem as well as the head waves, and the latter can be evaluated using the 
steepest descent method to provide a description of the tip diffracted body waves. It follows that all 
physically meaningful poles and branch points of $i(w + 9) must be located between the contours Ci 
and C2 and in this region the Sommerfeld amplitudes $i(w) should contain no physically meaningless 
singularities. Since inside the wedge we have \0\ < a, this means that all physically meaningful 
singularities — and only them must lie at a finite distance from the horizontal axis between the 
contours Ci and C2 as shifted horizontally by —a and a respectively, that is within the Malyuzhinets 
region 

{lo: ---a-2tan- 1 (e~ Im w ) < Re lo < — + a - 2 tan" 1 (e~ Im UJ )}. (2.8) 

To summarize, the scattered field can be fully described once an efficient algorithm is produced for 
calculating the Sommerfeld amplitudes $i(Lo) in the region (|2.8|) . We proceed with this task. 

3. Problem reformulation. In view of its symmetry with respect to the polar angle, the original 
problem naturally splits into "symmetric" and "antisymmetric" , corresponding respectively to the 
symmetric and antisymmetric parts of the incident wave. These involve functions ^>^(lo), such that 
we have 

$f(a;) = i[* i M + (-l) <+1 *«(-w)]. (3.1) 

We follow the approach pioneered in elastodynamics by Budaev (1995) and first substitute Ij2.6|l 
into the boundary conditions (|2.2|l and (|2.3(l to obtain the system of functional equations. We then 
employ the singular integral transforms to reformulate the problem as a system of algebraic and 
singular integral equations. 
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3.1. Functional equations. Let us start with the symmetric problem. The boundary conditions 
(I2T21) and (1231) imply 
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where a\(uj) — sin 2a;, (12(0;) — —cos 2w and 03(0;) = — 2cos 2 w. Introducing in the first terms of 
both (|3.2|) and (|3.3|) the new integration variable o>, such that cos w = 7 cos o>, dropping the check 
and transforming the contour of integration back to C the equations acquire the form 



f(uj)e [kr 



J duj = 0, 



(3.4) 



where, due to l|2.7|l . as | Im w\ — > 00, f(uj) — 0(exp(lm uj [2 — Re p])), with Re p > — 1. Using the 
Malyuzinets theorem (Malyuzhinets 1955) / is an odd trigonometric polynomial of the second order. 
It follows that the pair $+ (uj) satisfies ((3.2(1 and 1)3.3(1 if and only if it satisfies the functional equations 

*n{$o[ffM + »] + $ o"[ffM - "]} + t 12 {$+(cj + a) + $+(cj - a)} = Q+, 

t2i{*3"faf(w) + a] - $+[g(w) - a]} + t 22 ^+(uj + a) - $+(w - a)} = Q+ 

(3.5) 



where t\\ = cos 2w sin uj j \/j 2 — cos 2 w, £12 = £21 = sin 2u>, t 22 = —cos 2w; and we have 

Q+ = ct sin w + ct,sin 2w, 



(3.6) 



with ci, j', k = 1, 2 — unknown constants. The function 5(0;) = cos -1 (7 _1 cos w) relates the incidence 
shear angles to compressional reflection angles and its branch cuts are chosen so that the deformed 
contour of integration C = {uj = g{uj) '■ uj € C} may be transformed back to C without touching 
them. They are the segments 



[— 6h + 7m, Oh + Trn], 9^ = cos 7, n — integer. 

The branch of g(uj) is chosen so that it has the properties 

7T n 
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(3.7) 
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In order to investigate restrictions on c^ let us substitute expansions IB. 7(1 of the Sommerfeld 
amplitudes into ((3.5|) and equate the coefficients of the leading asymptotic terms in the resulting 
equations. Firstly, we note that in the symmetric case, l(B.7l) contains no constant terms and therefore, 
there are no exp(— 2\uj) terms in the left hand sides of these equations. This implies that cf 2 = c^ 2 = 0. 
Secondly, since the tip asymptotics of (j>f (kr, 8) contain the terms with the exponent 1 , these sides 
contain the exp(— iuj) terms. Equating the coefficients of the exp(— iuj) terms we obtain 
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Hence we have 



cJl = —dii tan a. 



(3.9) 
(3.10) 
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Note that if 7$^ + i$f m = 0, then c^ = -c£ = 0. 

It follows that the right hand sides in l|3.5|) might be - and as we show in §21 and are - nonzero, 
so that we have 

Q~l = c^sin lj, Q 2 = — cf tanasin u>, (3-11) 



where from now on, for simplicity, we use notation c[ = c^ . 

The antisymmetric problem can be treated similarly, with one minor modification: For all wedge 
angles a, expansion l)B.7(l might contain a nonzero constant term $q and therefore, the functional 
equations might contain a second order term. However, this term can be eliminated by subtracting 
<&^q from the Sommerfeld amplitude $jj"(w), redefining it in the process. The corresponding functional 
equations are 

t 21 {%[g(Lu) + a}+ 9q\0(w) - a]} + t 22 {$l(u + a) + $J" (w - a)} = Qf . 
tn[%[g(^) + a] -<f>o[g(uj) -a]|+ti 2 |$r( w + a ) - $1 {w ~ a)} = Q 



2 • 
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with 



Q x = c 1 sin to, Q 1 = c x tanasinw. (3.13) 



We note that the above reasoning involves only the asymptotic terms with p~ = or else with p^j = 1 
and N^ = 1, and therefore applies to all wedge angles under consideration. We note too that Budacv 
and Bogy have made several attempts to establish restrictions on the constants. They first used the 
arguments of the type outlined above in Budaev and Bogy (1998). By excluding from consideration 
the terms with p^ = 1, it is easy to reach the erroneous conclusion that all constants c fc are zero. 
In the static problems, such exclusion is justified, because the terms describe body translations. By 
contrast, in the dynamic problems, their presence is indicative of nontrivial phenomena. 

We proceed by discussing the singularities of the Sommerfeld amplitudes. First, we assume that 
the incident wave is plane or Rayleigh, so that it manifests itself in $ 4 (oj) in the form of terms which 
contain simple poles 9u in the strip | Re w\ < a. The functional equations H3.5(l and H3.12J1 can be 
recast as 



*j}=(ff(w) + a) \ =± ( ruM ria(w) \ / *|(?(w)-a) 
$f(w + a) / V r 2 i(w) r 22 (w) / \ $f(w-a) 



+ c? V7 2 -cos^ /e^)\ 
1 A(w) V e 2 (w) / 

where the reflection coefficients for the traction free elastic half space ijk(ui), j,k — 1,2 as well as 
the Rayleigh function A(lo) and functions ef (u>) are given in Appendix |jj The system (|3.14|) can be 
used to effect the analytical continuation from the strip | Re u)\ < a and thus find all poles 9u of 
the Sommerfeld amplitudes, with their respective residues, which are located in the strip Re u el — 
[it/2 — a, ir/2 + a]. The rationale behind the choice of the latter strip is clarified below. The poles 
are incidence and reflection angles of the respective incident, reflected and multiply reflected waves 
and their residues describe the amplitudes of these waves — see Budaev and Bogy (1995, Eqs. (17) and 
(18)). The first index in 9u refers to the mode of the wave and the second to its place in a sequence 
of all incident and (multiply) reflected waves (see Appendix [DJ) . 

Let us now again follow the above authors and introduce the decomposition 

$±H = $, ± ( W ) + $±(c), (3.15) 

where the unknown ^ (u>) is regular in the strip Re u> e I, and the known ^> t is 

$f(w)=^Res($f;0^y(w-0«), Re 9 a el. (3.16) 

t 

Above, an otherwise arbitrary function u(ui) should be chosen to be analytic everywhere inside the 
strip Re we I, except for a simple pole at zero, where it has the residue one. The weakest restriction we 
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can impose on behavior of <j{uj) at the imaginary infinity is that it grows slower than exp(| Im uj\k/2ol). 
Instead, we impose a stronger restriction — that it behaves as the amplitudes in (|2.7|l . If a(u>) possesses 
singularities which lie outside the strip Re we/ the functions <frf (uj) contain extra poles which describe 
waves that are outgoing at physical reflection angles but have nonphysical amplitudes. This causes 
no complication, since the corresponding singular terms in ^ i (uj) and &\ (uj) mutually cancel. 

Next we substitute decomposition (|3.15|) into the functional equations (|3.5[1 and then (|3.12|) to 
obtain the following inhomogeneous systems of equations for the regular components of the Sommerfeld 
amplitudes, 



[$+lg(oj) + a} + $+[g(uj)-a}}+B 

a{*Mu) +a}- $+[g(uj) - a]} + 
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and 
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where we use notations 



^21 (w) 

A= -^t-4 = -tan2w, 

t22{0J) 
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(3.19) 



To summarize, following Budaev and Bogy, the original problem can be reformulated as the 
following boundary value problem in the theory of analytic functions: Seek constants c x and functions 
®f {oj) , such that 

1. $j (w) are analytic for Re we/ and satisfy the asymptotic estimate (|2.7|) . 

2. the values that <&f (u>) take on the boundaries of the strip Re ujel are linked by equations 
(|3.17|) and l|3.18|l (that is, solve these equations for Re u — 7r/2). 

The above considerations and the properties of the Sommerfeld transform which are outlined 
in Appendix ^ show that such a pair exists if there exists a solution of the original problem. The 
uniqueness of $^ {uj) is a more complicated issue which we address in §5. 

3.2. Singular integral equations. Budaev (1995) has suggested to exploit the fact that for all 
functions F{oj) satisfying the first of the above assumptions, the singular integral transform 



(HF){oj) = ^-V.R 
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has the property 



H : F{u> + a) + F{lo - a) -* F{uj + a) - F{uj - a), 



Re W = ~. 



(3.21) 



This means that on the vertical line Re uj — n/2, the terms in the square brackets in l|3.17|l and l|3.18|l 
are related by H . The terms in the curly brackets are linked by a similar explicit transform, 



HF{uj) = ±VJ>. 



7r/2+ioo 



F{£)g'{m 



Rew = -, 



(3.22) 
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where (/(£) = dg/dt;. This suggests introducing new unknown functions 

X ± (w) = $f\g(w)+a]+$$\g{w)-a], Y ± (uj) = $f (w + a) + *f(w - a). (3.23) 

Note that the line Re u = ir/2 is of special significance, because the function g(u>) maps it onto itself. 
We can now use Ij3.21|l to transform Ij3.17|) and (|3.18(l into the system comprising algebraic equations 
and singular integral equations which hold on the vertical line Re u> = -f . This is the crux of Budaev 
and Bogy's approach. 

Changing to the new independent real variable n, such that uj = ir/2 + vq, the symmetric problem 
becomes 



where we have 



x + ( V ) + b( V )y + ( V ) = rt(v) - c t ^^^ L , (3-24) 

COSil Ti 

a(r,)Hx + (r)) + Hy + ( V ) = r+( v ) - c+ tima—^-, (3.25) 



x ± (r,)=X ± (-+ir 1 ), y ±( v )=Y ± {-+ir ] ). (3.26) 



Standardizing notations and substituting IJ3.24I) into (|3.25|) the problem transforms to a final singular 
integral equation in two unknowns, a function y + (rj) and a constant c{ , 

M+y + {ri) = q+(n) + 4qt(v), V - real, (3.27) 

where M + = TL — aTLb. Using the same approach the antisymmetric problem transforms to 

a(r])x-{r})+y-(r)) =rj(7?) + c i cosh2 > ( 3 - 28 ) 

M~x~(rj) = qo(rj) + c^q^(n), rj - real, (3.29) 

where M _ = 7i — Vria. The rest of the nomenclature can be found in Appendix [j] 

4. A New Numerical Schedule. Budaev and Bogy (1995, 1996, 2002) have advanced various 
implementations of the numerical schedule for computing $^(oj), all of which involve the following 
three steps: 

1. evaluating y + {rj) and x~{n) on the line n — (Re lu = 7r/2) by solving the singular integral 
equations (|3.27|l and (|3.29(l . and then evaluating x + (rj) and y~{r]) by solving the algebraic equations 
(TT2HI and d£2E); 

2. evaluating ^> i {to) in the strip Re we/, using the convolution type transforms H4.13(l and 
(|4.14|) below, with the kernels singular on the boundary of this strip; 

3. continuing the computed Sommerfeld amplitudes ^> i (w) analytically to the right of the strip 
Re uj e I by using the functional equations l|3.14|l . Recastin these equations to effect the continuation 
to the left of Re u> e, I. 

We have developed an alternative recipe for carrying out the first two steps. 

4.1. Solving singular integral equations in two unknowns on the line n = (Re u) = 

7r/2). Let us consider the symmetric case first. Operator M + is not analytically invertible, but Budaev 
and Bogy (1995) suggest that Eq. I|3.27|l can be rewritten as 

(Hd + K)y+(n)=q+( V ) + c+qi(v), (4-1) 

where 7i is the singular operator introduced above, analytically invertible in the space of bounded 
functions; K is a regular operator; and d(t) is an exponentially decreasing function. Importantly, Ti 
has the property 



oo 

fnf(r))dri = 0, (4.2) 
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and therefore, its range consists of all L 2 (R) functions, with the zero integral, where L 2 (R) is the 
space of all integrable functions of real variable. Budaev and Bogy (1996) state that they regularize 
(|4.1(l by applying Hr 1 to both its sides. They carry out numerical evaluation of the resulting singular 
integral equation by using (|4.2J) as a constraint and calculate c\ and y(rf) both at once. By contrast, 
below we argue that the right hand side of l|4.1|l belongs to the domain of TL^ 1 for only one value of 
C\ , carry out the regularization by finding this value and thus arriving at a singular integral equation 
in one unknown, 3/(77) . At present, our schedule works only for a < 7r/2. 

We start by observing that Eq. 1(4.1(1 is solvable only if its right hand side belongs to the range of 
TLd + K. We cannot describe this range explicitly. However, it is clear that {—Ky + + q^ + c^q^)(i]) 
should be in the range of TL. It follows that we must have 

00 

J [(Ky+)( V ) - q+( V ) - cfq+M] dr) = 0. (4.3) 

— OO 

All our numerical experiments confirm that neither gjj" (77) nor qf(r)) are in the range of 

TLd + K — by producing the nonzero "solution defects" Xq and A^~ defined by (|4.11l) below. Therefore, 
Eq. H4.1f> is solvable only if the right hand side of 1(4.1(1 is in the range. This gives the following 
relationship between cf and y + {rf): 



[{Ky+){rj)-q+{rj)]dr)[j q?(ri)d V \ . (4.4) 

— OO 

By substituting (|4.4|) into (|4.1fl . c± is eliminated and we obtain 

(Hd + P qt K)y+{ V ) = P qt q+, (4.5) 

where an unbounded projector 

00 00 

(P qt u)(ri) = Ufa) - q+( V ) J u(t) dt[J q+(t) dt] ~ ' (4.6) 

— 00 —00 

maps any function in L 2 (R) with a finite integral into the range of 7i and has the property 

ufa), for all u(r]), such that / u(t) dt = 0, 

(4-7) 

0, for u(rf) = q(rj). 

We refer to the function q{r[) as the projector kernel. 

Note that in l|4.5|l . the integrals of both P +Ky + (r]) and P q +qo(v) are zero, and therefore, the 
inverse operator 7i _1 can now be safely applied to both sides. Introducing on top of that a new 
unknown function y + (rj) = d 1 ' 2 {rf)y + (r}) the final regularized integral equation is 

y+( V )+L + y+(r 1 )=q + (r ] ), (4.8) 

where L + = d~ 1 ' 2 'H~ 1 P i . + Kd~ 1 ' 2 is an operator with a smooth kernel and the right-hand side q + (rf) — 
d~ ' (v)'^- 1 P„+ Qo (v)- The equation involves one unknown, y + (77) , and can be solved using a standard 
quadrature method (see e.g. Atkinson 1977). Note that normalizing the original unknown by d 1 ' 2 ^) 
rather than d(r]) leads to a new operator with a bounded kernel (c/. Budaev and Bogy 1995). The 
normalization achieves symmetrization of the kernel, so that whether 77 — > 00 or t — * 00, it exhibits 
the same singular behavior. 

Eqs. (|4.4|l and l|4.5|l imply 14.1(1 . This means that the combination of cf and a solution of l|4.5|l 
gives us the solution of (|4.1|l . However, in our code instead of solving (|4.8|l we implement a slightly 
different approach: Since qiiij) is rather complex, instead of P + we employ the projector P q2 , with 
the kernel 

©fa) = t r-7r-- ( 4 - 9 ) 
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Numerical experiments have shown that this kernel leads to a stable evaluation scheme. We then 
regularize and solve two equations 

(Hd + P q2 K)y+( V ) = P q2 qt(v), * = 0, 1 (4.10) 

(see Appendix IE)) . The "solution defects" 

oo 

A+ = f [(Ky+)(t) ~ qf(t)]dt, i = 0, 1 (4.11) 



turn out to be nonzero indicating that neither Qq(tJ) nor q{{r]) are in the range of Tid + K. It follows 
that the solution {y + {rj), cf) of Eq. (|3.27|1 can be obtained using 

4 = ~, y + (v) = vo(v) + 4vt(v)- (4.12) 

Ai 

The antisymmetric problem can be treated in a similar manner (see Appendix lEjl . 

4.2. Evaluating $^(w) in the strip Re en el. As already mentioned above, according to Bu- 
daev and Bogy's (1995), evaluation of $j (w) in the strip Re we/ can be carried out by using the 
singular convolution type integrals, 

oo 

*±( W ) = -L / /° {r,) dr,, (4.13) 

4m J cos ^ (£ - w + ixw)) 

witn /o + ('?) = 2/ + (»?) tanh 2r; - ic^cosh7y/cosh2ry + ir^ (ry)x'(ry), /^(r?) = ia: _ (»7)x'(?7); x'(??) = o?x/d?7 
and 

■DC 



^(")=^ / r ' , (4.U) 



■J, „. ^J_ /• /f(l)^ 

"4ai7 cos *L(?r_ w +iTj)' 

-oo z " z 

with fi(rj) = iy + {rj) and fi (r)) = a; _ (??) tanh2?7 — ic 1 "cosh?7/cosh2?7 — ir^??). 

If — as is the case for the Sommcrfeld amplitudes of the solution of the original problem — as 
Im co — > ±oo, the leading terms in l|B.7(l are 0(exp(±ipw)), with Re p > 0, then for a < ir, — Re p — 
ir/2a < and therefore, the above integrals converge. 

We start with the integral of the type l|4.14|l . Its generic form is 



oo 



— dr], | Re CI < a, (4.15) 



-OO 



cos ^(^ + i»?) 



where the new complex variable is £ = ir/2 — ui. When | Re £| < a, the integral l|4.15|l can be 
approximated using the trapezoidal rule. The approximation error is of order O [exp (— 2tt(t//i)], with 
h — the distance between the nodes of a uniform mesh and cr(u>) — the half width of the strip which is 
centred on the real line and inside which the integrand is regular. In our case, a < min{a — Re £, a + 
Re £} and as Re £ — > a, the accuracy of the trapezoidal rule deteriorates. Therefore, a more robust 
quadrature formula is required, with accuracy depending on the function f(rj) and not on parameter 
£. One such formula may be obtained with a modified sine function, 

2a ' 
Note that we have 

m»*) = { ; n ^ (4.i7) 
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and 

oo 

f(v)~ J2 f(nh)"h(v-nh) (4.18) 

n— — oc 

where the sum on the right interpolates f(r]). Therefore, the integral (|4.15|) may be approximated by 

oo 

i cos£t + ^ dT1 ^ Animnh) > (419) 

with the coefficients given by 

oo 

A n (0= I ^:f?^ . (4.20) 

J cos ^(f + 177) 

— oo 

These can be evaluated approximately by introducing new variables fj = r\ — nh and £ = £ + in/i. Then 
(|4.2Q(I can be rewritten as 



h [ sin \f] df) _ 

2a J sinh ^77 cos ^(£ + if/) 

OO 

e*^ 1 <2r/ /z /* e~^ J CZ77 



/7r7 /j f p — 7T 7 ? 1 *7ri 

(4.21) 



4 - Slnh ^a ^ C0S 5^ ^ + "^ 4 _ Slnh ^^ C ° S ^ ^ + "^ 

In the upper (lower) half plane where we can utilize the Jordan Lemma to evaluate the first (sec- 
ond) integral, each of the respective integrands possesses two sets of poles, zeros of sinh 7rf//2a, 
f) = ±2ami, m = 0,1,2 ... , with the respective residues 

± —p ( 4 - 22 ) 

7T COS ^£ 

and zeros of cos (tt(£ + if/) /(2a)), fj = i[£ ± a(2m + 1)], with the respective residues 

T— jr?— • (4-23) 

TV COS ^£ 

It is clear that a significant contribution to (|4.21|l is made only by the poles, 77 — and 77 = i[£ ± a] 
(in the upper and lower half plane respectively); other residues contain small exponential factors. 
Therefore, applying the Cauchy Residue Theorem and taking into account that the first pole lies on 
the contour of integration, we have 

3 ± ^* df] h (\ 



AaiJ^ sinh ^fj cos £ g + if j) C o S *-£\ 2 7 ' ^ j 

which — returning to the original variables £ and 77 — gives us a new quadrature formula 

00 

I Re £| <a. (4.25) 

Let us now consider the integrals of the type (|4.13|) . Its generic form is 

oo 

m 



I 



cos^K + ix(r/)] 



dr), | Re £| < a, (4.26) 
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where xiv) is a smooth monotone function. We could change the integration variable r\ to x(v)i reduce 
(|4.26(l to the integral of type (|4.15|l and evaluate the result using a uniform mesh in x- However, both 
integrals (|4.15|l and (|4.26(l involve the solution of (|3.27|) . and therefore, it is more reasonable to 
evaluate both integrals using the same mesh. Then following the same reasoning as above, (|4.2f>fl may 
be approximated by 



I WITT- ( ^ ra £ ®n(Of(nh), (4.27) 

J cos ^-[^ + i X (ry)] ^ 

where the coefficients are given by 

oc 

J COB ^k^ + 'xiv)} 

and the main contributions to (|4.28|l are made by the zero 77 = (77 = nh) of the hyperbolic sine in 
LOh and the zero f\ = \a± — nh (n = ia±) of the cosine- function, where £ + ix(io±) — -foe- The latter 
equation implies that ix(ia±) = — sin - (7 _1 sin a±) = =Fa — £, and therefore, we have 

a± = sin" 1 [7 sin (£ ± a)] , (4.29) 

with Re a + > and Re a_ < 0. Applying to (|4.28l) the Cauchy Residue Theorem and noting that 
the pole f\ — lies on the contour of integration we obtain 



oc 



h f e ± ^ 1 dfj 



4ai J sinh ^fj cos ^ [£ + i X (f) + nh)} 

r 1 e Ti(a±+'mh) . 

/l l ± 2cos £[£ + i X (nh)} T X '(a±)sin £(a ± + inh)i' {4M) 



where x'( a ±) — cos a±/\/ / j 2 — sin a±. Returning to the original variable 77, the resulting quadrature 
formula is 



00 

/■£ I \ ! "^ 1 

cos « 'Ufo)]* 7 " \§Jc^«: 






K + ix(nfc)] 



(-i)"[ r? , - e I, M - M + r? r e ^ ^- d }^ nh ^ 

X (a+)sm ^(a+ + mft) x'(a_)sm ^(a_ + m/i) J J 

| Re C| <a. (4.31) 

The first terms on the right of both l|4.25[) and i|4.31|) effect the trapezoidal rule and the second give 
a correction. 

5. Code Testing. Using the above considerations we have developed a new code for evaluating 
the Rayleigh reflection and transmission coefficients for elastic wedges (see Appcndix[UJ). The integral 
equations we solve have the form of the Frcdholm equations of the second kind, but it can be shown 
that the operators involved are not Fredholm (c/. the statements in Budaev and Bogy 1995, p. 251). 
We possess no analytical proof that these equations can be solved uniquely. Nevertheless, our code 
produces a solution, and below we describe verification tests that allow us to state with confidence that 
when transformed back to the physical space this solution satisfies the original diffraction problem. 
We also describe successful validation tests, comparing output of our code with published numerical 
and experimental data. Of course, the positive outcomes of these tests do not constitute a theoretical 
proof that the code is correct. Note that throughout this section we characterize materials by their 
Poisson's ratio v, where 7 = y^I — 2i/)/[2(l — v)}. 

5.1. Code verification. We have designed verification tests to establish that the computed 
functions ^i(co) are the solutions of the original physical problem, in particular, that they 
(i) are bounded at imaginary infinity; 
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(ii) are analytic at the boundary of the strip Re uj el; 
(iii) possess only physically meaningful singularities. 

The property (i) is confirmed by direct examination of the computed functions x~ (rj) and y + (77) 
divided by exp(— \r]\). At large \-q\ the ratios appear to behave as 0(1). It follows that the amplitudes 
$i(w) obtained by the analytical continuation must be bounded at the imaginary infinity, rj = Im to — > 

~Xj. 

Since the last step in the analytical continuation is carried out strip by strip, all 2a wide, there 
is no guaranty that any computed &i(uj) should be smooth at the boundaries of the initial strip 
Re to el. However, examination of the numerical output related to Figs. 15 . II - 15 . 3l confirms that our 
approximations are smooth: The numerical derivatives of computed <&i(uj) jump at the boundaries of 
the strip by about 10 . It follows that the property (ii) is satisfied. Interestingly, when attempting to 
solve an incorrect problem, with the constants c* put to zero, the computed ^(w) themselves jump 
at the boundaries by about 10~ 2 . 

Similarly to (ii), the property (iii) should be satisfied by the Sommerfeld amplitudes of the so- 
lutions of the original wedge diffraction problem, but it is not obvious that the computed solutions 
of the corresponding functional equations should satisfy it as well. Indeed, the way they are con- 
structed assures that the computed amplitudes $i(w) have physically meaningful poles and since the 
functional equations that are used to effect the analytical continuation involve reflection coefficients 
ijkioj), with the branch point at u> = 9] x and poles at ui — ±i/3ri (sec (I3.14|) and Appendix 0J , they 
possess physically meaningful branch points and Rayleigh poles too. However, by the same token, 
the analytical continuation scheme might endow these amplitudes with extra, physically meaningless 
singularities. Remarkably, when the incident wave is compressional or shear at lo = — a ± i/3jj, all 
our computed residues are of order 10 -7 , i.e. are numerical zeros. Thus, the computed Sommerfeld 
amplitudes possess no Rayleigh poles corresponding to physically meaningless Rayleigh waves incom- 
ing from infinity. Furthermore, Fig. I5.1H5.2I which respectively relate to a purely symmetric and 
purely asymmetric case, confirm that inside a neighborhood of zero which includes the strip Re to el, 
the computed Sommerfeld amplitudes possess the symmetries described in (|3.1|l . that is <&q{lu) and 
<&i(uj) are odd while $>*((*}) and $^~ (w) are even. (Outside this neighborhood, the symmetries are not 
apparent in Figs. l5"lTl5.3l due to the accumulation of numerical errors.) As we show in Appendix IFI 
such symmetries imply the absence of physically meaningless branch points. 

(a) (b) 



___ U=J _ =- ; 




V 



Figure 5.1. The computed Sommerfeld amplitudes: (a) Re <3>o(^) — dashed line and Im &o{^) — solid line, (b) 
Re "3>i(w) — dashed line and Im <I?i(aj) — solid line. Wedge angle 2a = 70°, / = [0.96,2.18], Poisson's ratio v = 0.25, 
incident wave — compressional and 9 mc = 0° . 



The properties (i) , (ii) and (iii) of the computed Sommerfeld amplitudes respectively imply that 
they possess all the properties expected of the Sommerfeld amplitudes of the solutions (pi (kr, 9) of the 
original diffraction problem, so that their corresponding Sommerfeld integrals satisfy (i) the Helmholtz 
equations and correct tip condition; (ii) zero stress boundary conditions; and (iii) radiation conditions 
(which exclude nonphysical Rayleigh or head waves incoming from infinity). 

Fig. 15. II provides one more confirmation that the computed functions "^(cj) are the Sommerfeld 
amplitudes of solutions 4>i (kr, 9) of the original wedge problem: They show that for the symmetric com- 
pressional wave incidence, both $i(ui) are imaginary and therefore, the corresponding displacements 
are real. This is consistent with the physics of the problem, since unlike with the symmetric shear 
wave incidence, there is no total internal reflection, that is no imaginary displacement component. 
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(b) 
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Figure 5.2. The computed Sommerfeld amplitudes: (a) Re ^(w) — dashed line and Im "J>o( w ) — so/id ime, ffcj 
Re <&i(loi) — dashed line and Im "J>i(u;) — soKd line. Wedge angle 2a = 70°, I = [0.96,2.18], Poisson's ratio v = 0.25, 
incident wave — shear and 9 mc = 0° . 

(a) (b) 




Figure 5.3. The computed Sommerfeld amplitudes: (a) Re "^o^) — dashed line and Im Qoi^) — solid line, (b) 
Re Q\(u)) — dashed line and Im <&\(lo) — solid line. Wedge angle 2a = 70°, I = [0.96,2.18], Poisson's ratio u = 0.25, 
incident wave — Rayleigh and Re 8' nc = —a. 



Finally, a numerical stability of the scheme is ascertained by the fact that different choices of 
adjustable function u(u)) in IJ3.16JI all give similar results (see Appendix |HJ). 

5.2. Code validation. Our first validation results are presented in Tables I5TTI and ET21 where the 
approximate values of amplitudes and phases of reflection and transmission coefficients R and R tian 
as computed with our code are compared with numerical results of Fujii (1994). Each Fujii's column 
contains values corresponding to different choices of an adjustable parameter. The parameter allows 
one to evaluate singular integrals on the real axis by moving the poles away from the axis into the 
complex plane. This is equivalent to employing the radiation condition at infinity in the form of the 
limiting absorption principle. From the physical point of view, the singularities cannot be moved too 
far. However, when they are too close the evaluation algorithm becomes numerically unstable. The 
top rows in the tables are obtained with the parameter values that correspond to a more physically 
meaningful situation and the bottom ones, with the values that give better numerical stability. The 
tables demonstrate that for the larger wedge angles the agreement with our computations is quite 
good, but for the smaller ones our values lie outside Fujii's range. This is not surprising, because 
when the wedge angles are small there are many multiply reflected waves and the residues of many 
resulting poles are large. For this reason, when the wedge angles are small, the present version of our 
code looses its numerical stability. 

To continue, in Figs. 15.41 (a) and (b) we present our Rayleigh reflection and transmission co- 
efficients as functions of the wedge angle, computed for v = 0.234. They fit Fujii's numerical and 
experimental data extremely well (see our Fig. 15.51 or Fujii 1994, Fig. 7). Note that on Fujii's plots 
the solid lines represent his numerical results and discrete points, his experimental data.) Indeed, for 
the wedge angles between 45° and 150° we cannot put the results on the same graph — there is no 
visible difference. Note that the jumps in the phase of the reflection coefficient that take place at 
the wedge angles of about 45° and 145° are from 180° to —180° and —180° to 180° respectively, and 
therefore, no jumps in physical quantities take place. For the wedge angles between 150° and 180° 
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the reflection coefficients are practically zero. This is understandable, because when the wedge angle 
is 180° there is no reflection. In this region, the phases of our reflection coefficients differ from Fujii's 
but the limiting value of 90° agrees with the one obtained by Gautesen (2002a). It appears that in 
this region Fujii's scheme looses its stability. 

Table 5.1 
Rayleigh reflection coefficients computed with our code and Fujii' (see Fujii 1994)- v = 0.25. 



wedge 




|i? ref | 


arg 


R re{ 


angle 


Fujii 


This paper 


Fujii 


This paper 




0.50552 




-169.87° 




50° 


0.49924 




-169.54° 






0.49278 


0.47427 


-169.07° 


-161.4° 



150° 



0.05257 
0.05252 
0.05236 
0.05217 
0.05196 
0.05151 



0.05197 



170.53° 
170.14° 
169.85° 
169.65° 
169.53° 
169.43° 



170.5° 



Table 5.2 
Rayleigh transmission coefficients computed with our code and Fujii' (see Fujii 1 



v = 0.25. 



wedge 


| 


h" 


ran 1 


arg 


Dtran 


angle 


Fujii 




This paper 


Fujii 


This paper 




0.49123 






-32.59° 




50° 


0.48372 






-32.84° 






0.47552 




0.55189 


-33.00° 


-26.9° 



150° 



0.78940 
0.78899 
0.78866 
0.78842 
0.78825 
0.78800 



0.78942 



52.79° 
52.80° 
52.81° 
52.82° 
52.84° 
52.86° 



52.9° 



The amplitude curves reported by Budaev and Bogy (2001) are the same as ours (see ibid, Fig. 
6 and our Fig. 15.41 (a)), but for the larger wedge angles, the phase of their reflection coefficient is 
somewhat different — see Fig. 15.41 (b). The discrepancy might not be crucial, because at these angles 
the amplitudes of the reflection coefficients are very small, but the problem is indicative of numerical 
instability. Note that the results on the wedge angles greater than 180° as presented by Budaev and 
Bogy (1996) are incorrect — see their errata (Budaev and Bogy 2002). Note too that even though 
Poisson's ratio used by Budaev and Bogy (2001) is v = 0.294 the above comparison is valid: The 
coefficients should not be effected by a small difference in v (see e.g. Fig. I5.6h 

We finish this section by comparing our computed Rayleigh reflection and transmission coefficients 
for the quarter space with Gautesen's. On taking into account that Gautesen's coefficients are complex 
conjugates of ours and thus, our phases must have the opposite sign, the agreement between the 
calculations is very good (see Fie. 15.611 . 

6. Conclusions. We have studied the properties of the underlying integral operators and devel- 
oped a new numerical schedule for solving the singular integral problem that arises in Budaev and 
Bogy's approach to diffraction by two dimensional traction free isotropic elastic wedges. We have also 
developed new quadrature formulas for evaluating the singular convolution type integrals that are 
utilized in this approach. Although the analytical justification of the method is not entirely rigorous, 
the code has undergone a series of stringent internal verification tests directed at establishing that 
it solves the original physical problem as well as validation tests against numerical and experimental 
results reported by other authors. It appears to be successful when simulating diffraction by wedges 



16 

(a) 



. KAMOTSKI et al. 



(b) 




60 u 



180 u 



120 u 
WEDGE ANGLE WEDGE ANGLE 

Figure 5.4. Rayleigh transmission coefficients (solid line) and reflection coefficients (dashed line) computed with 
our code versus the coefficients computed with Budaev and Bogy's code (squares and circles respectively — see Budaev 
and Bogy 2001, Fig. 6). Poisson's ratio v = 0.234, incident wave — Rayleigh and Re 9 mc = —a. 

(a) (b) 
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180* 



Figure 5.5. 



The Fujii's computed (solid line) and experimental (dots) transmission and reflection coefficients: 
(a) amplitudes, (b) phases. Poisson's ratio u = 0.234, incident wave - Rayleigh and Re mc = —a. Reproduced from 
Fig. 7 in Fujii (1994). 

(4 , 0*°, 
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Figure 5.6. Rayleigh transmission coefficients (solid line and squares) and reflection coefficients (dashed line and 
circles) computed respectively with our code and Gautesen's code (see Gautesen 2002a, Figs. 3 and 4): (a) amplitudes, 
(b) phases. Wedge angle 90°, incident wave — Rayleigh and Re mc = —a. 



of angles between 40° and 178° of plane incident waves, compressional or shear. When the incident 
wave is a Rayleigh the lower limit of applicability goes up to 45°. 
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Appendix A. Analytic properties and asymptotic expansions of $i. 

In this Appendix we first follow Malyuzhinets (1957), Osipov and Norris (1999) and Budaev and 
Bogy (1995) and explain in more detail why the elastic potentials fa may be represented in the form 
of the Sommerfeld integrals. We then elucidate the analytic properties and asymptotic expansions of 
functions 3>j(w). 

A.l. The representation of fa in the form of the Sommerfeld integral. We develop our 
arguments for the potential fa only. The potential fa may be treated in a similar manner. Let us 
introduce the Laplace transform 



*i(«,0) = fc/ ip{kr, 9)e~ ksr dr. 
o 



(A.l) 



The behaviour of il>(kr, 9) at infinity is known from the physics of the problem: The potential comprises 
the incident and reflected waves, waves diffracted by the wedge tip, head and Rayleigh waves. As 
r — > oo the corresponding terms in the asymptotic expansion of the function ip are of order 0(1). 
Their derivatives in both r and 9 are of the same order. 

It follows that the function <I?i(s,#) is defined and regular on the whole half-plane Re s > 0. As 
s — ► vr + 0, where r is real, the function $i becomes the Fourier transform of ip{kr, 9). For any 9, the 
function $i(ir + 0, 9) has a finite number of singularities. Each is associated with one of the waves 
mentioned above. Excluding these singularities, the function $i(ir + 0, 9) is an analytic function of r. 

Using the inverse Laplace transform, ip(kr,9) may be written as 



ip(kr,0) 



1 
2tti 



cr+ioo 



*i(*, 



^krs 



ds. 



(A.2) 



Introducing a new independent variable to, such that s 
integral 



i cos u>, Eq. (|A.2|I becomes the Sommerfeld 



il>(kr,0) = I $i(w,60e L 
c 



krcos 



'duj, 



(A.3) 



where we use the notation 



$i(w,( 



2tt 



sin lo <J?i (i cos to, 9), 



(A.4) 



and C is a (J shaped Sommerfeld contour running from ioo to it + ioo (see Fig. IA.1|) . 



Figure A.l. The Sommerfeld contour. 



Since <l>i(s,#) is regular in the right half-plane Re s > and the new variable s — icos to maps 
this half-plane onto the half-strip {0 < Re lo < n, Im lo > 0}, it follows that $i(a>, 9) is regular inside 
this half-strip. 
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A. 2. The analytic properties of $1. The tip conditions assure that as kr — > 0, ip(kr,8) and 

derivatives dgip, and dgip remain bounded. The boundedness of ip implies that for a fixed args, such 
that — 7r/2 < args < n/2, the function $i(s, 9) — 0(s^ 1 ). Taking into account the formula (|A.4|I this 
in its turn implies that 4>i(o;, 9) = 0(\). The same estimate is valid for the derivatives of <3E>i, so that 
deK $*i = 0(l). 
Let us assume that 

1. excluding a finite number of poles and branch points on the boundary of the half-strip {0 < 
Re lo < it, Im lo > 0}, for any fixed 9, the function $i is regular inside the half-strip 

— e < Re id < e + ir, Imw> — e; (A. 5) 

2. $! is regular in the vicinity of 0; 

3. for large Im lo, inside the half-strip ljA.5|) . <f>i, dg&i, dg9& 1 = 0(1). 

These assumptions may be justified using results of Kamotski and Lebeau (2006); they allow us to 
transform the contour C (see Fig. lA.lfl into the contour 

{Re lo — —e', Im w> const, 

-e' < Re u < tt + e', Imw= const, (A.6) 

Re w = n + e', Im lo > const, 

where the constant is sufficiently large and < e' < e. 

Following Budaev (1995) let us transform the contour C in the integral (|A.3J) into the contour C e > 
and substitute the resulting integral into the Hclnrholtz equation. We obtain 

/ (pi -dfj^^Lo, 9)e kircos u dw = 0. (A.7) 

The Nullification Theorem implies that inside the contour C e > we have 



(dt- %)&i(w,e) = o. (a. 



Excluding the poles and branch points, the function $i and its derivatives are regular functions inside 
the contour l|A.6|) . Therefore, (|A.8J> is valid inside the region l|A.5(l . 
Eq. (|A.8fl implies 

$i(w,0) = ^i 1 ^^ + 9) + ^\lo - 9). (A.9) 

Let us apply the operator d^ + dg to both sides of (|A.9|I . This gives 

(d u + dg) ^(lo, 9) = 2[<S>[ 1) }'(lo + 9). (A.10) 

Excluding its poles and branch points, the function in the left-hand side of IjA.lOfl is regular inside 
the region {|0| < a, Im lo > —e, — e < Re lo < tt + e}. Therefore, both $^ ' and [<&\ ']' are regular 
inside this region and excluding its poles and branch points, &\ (lo) is regular in the half-strip 

— a — e < Re lo < e + a + tt, Im lo > —e. (A. 11) 

Applying the operator d u — dg to both sides of (|A.9|I we obtain 

(flu - d 9 ) fci(w, 9) = 2[<S>[ 2) }'(lo - 9). (A.12) 

Eq. (|A.12J| implies regularity of $\ (lo) in the half-strip 

— a — e — 7r<Re lo < e + a, Im lo < e. (A. 13) 

As kr — » oo, regularity of <&\ and $^ allows us to evaluate the integral (|A.3J| using the stationary 
phase method. The stationary points are lo = and lo = tt. The asymptotic term $i(0, #)exp [i(kr + 
7r/4)]/v27rfcr, evaluated at lo = is zero, because there can be no cylindrical wave incoming from 
infinity. Therefore, we have 

$i(O,0) - *P(0) + $i 2) (-6») - 0, (A.14) 
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and in the neighbourhood of 9 = 0, 

It follows that (|A.15|I is satisfied identically. Therefore, we have 



,(D 



(!)/ 



#i(«, 9) = $Y (w + 0) - <5>\>{-lj + 



and 



iP{kr, 9) = / [*i(w + 0) - *i(-w + 0)]e fcircos " dw, 
where we dropped the superscript (1). A further change in notations finally gives us Eq. 1)2. 5[l . 



19 

(A.15) 
(A.16) 
(A.17) 




Figure A. 2. The regularity domain o/<E>i(u>). 



A. 3. The domain of regularity and location of singularities of $i(w) . Let us return to 
(IA.10|I . Its left-hand side is regular in the region {Im u> > 0, < Re ui < 7r}; \9\ < a}. Therefore, 
the function < E , i(w) has no singularities in the half-strip {—a < Re ui < it + a, Im ui > 0}. It can be 
similarly argued that the half-strip {— n — a < Re u> < a, Im lo < 0} contains no singularities either 
(see Fig. IA.2I where the dashed and solid lines are the boundaries of the two half-strips in which 
the function $i(w) is holomorphic.) However, there might be singularities on the boundaries of these 
half-strips. 

Let us consider the solid lines 

{Re lo = -a, Im w > 0} U {Im w = 0, -a < w < +a, } U {Re u> = a, Im u> < 0}. (A. 18) 

The singularities which lie on these lines give rise to incoming surface waves. Let us discuss this point 
in more detail. Let us consider the Sommerfeld integral 



<J>i(w 



^l/crcos u 



duo. 



(A.19) 



cue 



When evaluating (|2.5(l . any pole ujq of $i(w), such that —a < loq < a, gives rise to a plane wave with 
the phase factor exp [ifcrcos (wo ~ 0)\ (see (IA.20|I ^. that is a plane body wave propagating along the 
ray 6 — luq which lies inside the wedge. By the same token, any pole uo\ = —a + i/3^ of <i>i(w), with 
fin > gives rise to a plane wave with the phase factor 



^ifcrcos (ui—9) 



^ifcrcos (a+0) cosh/^^p— /crsin (a-\-9) sinh f3n 



(A.20) 



so that its amplitude is exponentially small everywhere except for a small neighbourhood of the wedge 
face 9 = —a. Thus, the pole io\ = —a + i/3r, (3r > gives rise to the surface (known as Rayleigh) 
wave which propagates from infinity to the wedge tip along the wedge face 9 = —a. Similarly, any 
pole 0J2 — ol — i/3fi corresponds to the Rayleigh surface wave which propagates from infinity to the 
wedge tip along the wedge face 9 = a. 
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A. 4. The asymptotics of $^ at infinity. Let us finish this Appendix by discussing the be- 
haviour of < &i(cj) at infinity. Again we carry out the argument for $1 = <&i only. The Laplace 
transform $i(s, 0) has an asymptotic expansion similar to l|B.2|) . with the expansion coefficients that 
are functions of (Fedoryuk 1977 and and Osipov and Norris 1999, Theorem 2). Significantly the 
p = terms are either absent or involve constant coefficients (Kozlov et al. 2001). Therefore, the 
function 



* / 1 (w) = i[a ta) * 1 («- 



■$) + 9e§i(w + 6)] 



)=0 



(A.21) 



has the expansion of the type (|B.2|) without a constant term. Using l|A.4(l this implies that as Im to — ► 
co, ^i(w) — > and therefore, $i(u;) — > const. When the constant is chosen to be zero (or another 
specified value) we call the corresponding function $i(o>) the Sommerfeld amplitude. 

Appendix B. The tip asymptotics of the elastic potentials and asymptotics of the 
Sommerfeld amplitudes at infinity. The behavior of solutions of the elliptic problems in regions 
with piecewise smooth boundaries has been studied by many authors (see Nazarov and Plamenevskij 
1994 and references therein). A rigorous theory has been developed after a breakthrough by Kon- 
drat'ev (1963) who constructed and justified the field asymptotics in the vicinity of conical and edge 
points. The theory implies that the solution of the underlying Lame problem must have the asymptotic 
expansion 



u(kr, 9) 



l,m=0 



{kr) q 



N m -1 
n=0 



Uf. 



,(0) (Infer)™, fcr-^0, 



(B.l) 



where for one m, q m — and N m = 1 (otherwise, the tip conditions are violated); and for any other 
to, Re q m > and q m is a root of a transcendental equation, with a natural number N m being its 
multiplicity. 

The asymptotic expansions l|B.l|l can be differentiated and substituted into the boundary con- 
ditions. Therefore, using (|2.4|) similar expansions may be written for the elastodynamic potentials 
t/>f(fcr,0)as 



N= 



ct>t(kr,6)~ £(fcr)»*+ M £ 



i.l.m.n 



(0) (Infer)™, kr^O. 



(B.2) 



e,m=o 



n=0 



Let us arrange the sets of exponents {p+ , m = 0, 1, ...} and {p n 



0, 1, ...}, each in order of 



the increasing real part. Following Kozlov et al. (2001), these sets may be described as follows: Each 
contains 1, while any other element is a root of the transcendental equation 



and satisfies condition 



(p ± + 1) sin 2a ± sin 2a(p ± + 1) 



Re p ± > -1 







(B.3) 



(B.4) 



(otherwise, the tip conditions are violated). Note that if in l|B.l|) . a q m =^ 0, then in (|B.2|) the 
corresponding p+ or j?~ equals q m — 1, but applying the nabla operator to the term with q m = and 
£ = always gives us zero — because the corresponding N m = 1. Thus, for q m = 0, only the next, r 2 , 
term in (|B.1|I gives rise to a nonzero term in (|B.2|) . The corresponding p^ = 1 . Note that while for 
any wedge angle there exists an to, such that p~ n = is a solution of (JB.3J1 . in our range of wedge 
angles 2a £ (0, 7r), all p+ differ from zero. The full set of solutions of the transcendental equations 
HB.3I) is described in Kozlov et al. (2001). The main facts can be summarized in plots representing 
roots of the transcendental equations as functions of the wedge angle 2a (see e.g. Ting, 1984). The 
analysis of these plots shows that in our range of wedge angles, the tip conditions are assured for 
those non-negative exponents with the minimal real part that are either 1 or else are solutions of the 
corresponding transcendental equations, with the real part less or equal than 1. In other words, all 
leading exponents in l|B.l|l . that is, the exponents with the minimal real part lie in the strip 



< Re p„ < 1. 



(B.5) 
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The root loci in Ting (1984) also show that at one wedge angle, 2a» w 0.87T, we have a degeneracy: 
In the corresponding symmetric problem, the exponent with the minimal real part, p+ « 0.76 is 
a multiple root of the corresponding transcendental equation IJB..3J1 . The corresponding multiplicity 
N+ = 2. There are no multiple roots p$ which have the minimal real part and simultaneously satisfy 
(IB.5P . It follows that for a* , the leading terms in l|B.l|) are 



(f>t 011 (kr) p i\iikr = O((kr) p ), < p < p+ . (B.6) 

The behavior of potentials (f>\ (kr, 9) in the vicinity of the wedge tip dictates the asymptotic 
behavior of the Sommerfeld amplitudes 3>f (ui) at infinity: For example, it is easy to check that for 
any small e > 0, as Im u> — > oo, ^^(w) have expansions 

*f (w) ~ J2 ^^ + 0(e i(1+£) -), « ^ a„ < a < J, 

0<Ro p±<l 

$+( w )^$»we ip ^ + $ l *e ip ^ + $ l ie iw + 0(e i ( 2+e ^), a = a,. (B.7) 

Note that in the symmetric case, the expansions contain no constant terms (so that the leading 
exponents are 1 and possibly a solution of (|B.3J) . with the real part in (0,1)), but these may be 
present in the antisymmetric case (so that the leading exponents there are and 1). 

Appendix C. Radiation conditions at infinity. 

In wedge problems, the radiation conditions at infinity are readily expressed in the form of the 
limiting absorbtion principle: As kr — > oo, the diffracted waves u dlff with Im k < should decay 
exponentially. The solution of the original problem is obtained by taking the limit as Im k — > 0. 

Recently a version of the limiting absorbtion principle for the elastic wedge has been justified in 
Kamotski and Lebeau (2006). As a result, the new radiation conditions have been introduced which are 
consistent with the physics of the dynamic problem under consideration. The new radiation conditions 
have allowed the authors to prove that the elastic wedge problem stated in terms of displacements has 
a unique solution. Reverting to the language of potentials the Kamotski and Lebeau conditions may 
be re-written as 
the inner radiation conditions 

\i-Ykcf>£ S + d r <t)f e \ 2 rde^O, (C.l) 

\ikcpf s + dr^l 2 rdO^O, fcr^oo (C.2) 



\<a—(kr)~ 1 + s 



\<a~(kr)~ 1 + s 

surface radiation conditions 



/ 



\k R 4>f n + d r 4>f n \ rd<9^0, r^oo, (C.3) 



a-(kr)- 1 + '<\0\<a 



where ka = Q/cr plus the requirement that as kr — » oo the diffracted field satisfies the boundedness 
condition. Following the methods of Kamotski and Lebeau (2006) it must be possible to show that 
the same conditions are sufficient for proving uniqueness of the elastic wedge diffraction problem stated 
in terms of elastodynamic potentials. 

Appendix D. Geometrico-elastodynamic poles of $;. 

In this Appendix we present recursive formulae for evaluating the geometrico-elastodynamic poles 
of $j(w) with the corresponding residues. 

D.l. The symmetric problem. Given an incident P (compressional) wave 

$ lc (kr, 9) = e ^ krcos ( e - e o nc )l ; ( D .l) 

the symmetric function <&q (cu) has two poles at the angles of incidence #oi and #02, with the corre- 
sponding residues Res $J(#c«)) ^ = ^'^' These are 

n nine n flinc 

Res$+(M = -^r- (D.2) 
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The poles lie in the strip | Re ui\ < a. It follows from the first equation in (|3.14|l that $o"(w) also has 
P-P poles Bog, £ = 3,4, with the respective residues Res $^(0of)- These describe the once reflected 
compressional wave and are given by 

9oi = &0J-2 + 2a, 

Res $+(#,«) = r n (g- l (eo,t- 2 + a))Res $J(0 o ,*-2). (D.3) 

The S (shear) waves generated on the first reflection are associated with the P-S poles to — 9u, £ = 3, 4, 
with the corresponding residues Res $J~ (^i^). These are given by 



ni 



-i, 



h,e-2 + a) + a 



, - ,u. , ._ r 2i(g H^-2 + a)) „, , Ti - , 



Res S+(0 U ) = — ^ ~ , - Res *J(^-a). (D.4) 

3 (t/w - a) 



Similarly, the poles 9 i and On, £ — 5, 6 of the type P-P-P and P-P-S are described by formulae ID. 3(1 
and ID. 41) respectively. The P-S-P poles 9 g, £ = 7,8 are given by 

9 M = g{6i,t-A + a) + a 

Res ^(6 i) = ri 2 (0i / _ 4 + a)</(0 M _ 4 + a) Res §i(0 M _ 4 ), (D.5) 

and P-S-S poles U ,£=7,8 by 

9u = 9i,£-4 + 2a, 

Res $+(#«) = r 22 (0 M _ 4 + a) Res $+(0 M _ 4 ). (D.6) 

The above formulae should be applied recursively until one evaluates all poles within l|2.8|l . Due 
to the symmetry of $>q(uj) and $|"(w), — 6 are their respective poles as well, with the corresponding 
residues Res $J(— 9 e) = Res $d(^) and Res <&~l{—9u) — —Res ^{9u). If the recursive procedure 
is initiated outside the strip | Re uj\ < a, then, as equations 1)3.14(1 suggest, other pairs of poles appear 
to be possible. However, direct calculations show that such points give a zero contribution. 

D.2. The antisymmetric problem. The antisymmetric problem allows a similar treatment, 
the analytic continuation should also be performed using Eq. ((3.14|) with the "— " superscipt. 

Appendix E. The integral equations for one unknown. 

Symmetric problem. The two final regularized integral equations to solve are 

y+( V ) + L + y+(v) = q+iv), * = 0, 1, (E.l) 

where L + is an operator with a smooth kernel 



1 f l + (r],t)t&nh2t 

/ W)\/d(ji) cosh ^r] 



(L+u)( V )=-— J 7 =^=f= — u(t)dt, (E.2) 



with 

cosh^-r r tanh2£ tanh2r 



r cosh tj-t ( 

J+(„,t)=V.R / . ^ { — 

J smhJ-{T-rj) l X ( 



^ sinh^(r - 77) l X '(t)sinh^(i - r) s inh^[ X (i) - x(r)} 
The respective right hand sides of Eqs. ((E.l(l are given by 



dr. 



1 7T 

1 r coshT^-i 

qtM = = — V.R / .,*%* , Qt(t) dt, (E.3) 

2a ~°° 

with q[{rf) given in Appendix 0J On solving l(E.l(l . y + {rf) is obtained using 

y+{r,) = d-^ 2 ( V )[y+( V ) + c+y+( V )}, (E.4) 
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where c^ = — \q /X^ and we have 

[yi{t)d- l l 2 {t){Bn){t)-qf{t)]dt, * = 0,1, 

OO 

7 2 tanh2t f sinh2(x(i) + r) dr 



(B+I)(i) 



2ai 7 l + 27 2 sinh 2 (r + x(i)) sinh 27i T 



Antisymmetric problem . Analogously to the symmetric case, the two final regularized integral 
equations to solve are 

xT{r,)+L-xT{r 1 )=qr{m), i = 0, 1, (E.5) 

where the integral operator is 

OO 

(L-u)( V )=-— / )A4 «t)dt, E.6 

4a 2 J / o ^d(iyyd(^)cosh^x(»7) 

with 

cosh 2a" r / tanh2t | tanh2x~ 1 ( T )(x~ 1 )'(T) \ , 



sinhXfvM-rli^hX 



^sinh^W^-rllsinh^W^-r) S inh^[ X -i(r) - 1] 
The respective right hand sides of Eqs. (|E.5|I are given by 



oo _ 

i r cosh T^-t 

2aiVd(^)cosh^-x('7) •/ swh &-XW] 

with q~ (rj) given in Appendix |jj As above, x~ is obtained on solving l|E.5() using 

x~(rj) = d- l / 2 (T!)[xo(v) + ^x-(n)}, (E.8) 

where c^ = — A^/A^ and we have 

[xr{t)d- l l\t){B-^){t) - q r(t) x '(t)]dt, i = 0,1, 

, _ .... tanh2t f 00 tanh2(£ + r) dr 
(-« X )(*) ~ 



2ai 7_oo sinh ^r 

Appendix F. The branch points of $j(w). 

Let us show that the Sommerfeld amplitudes <&j(u>) that satisfy the functional equations (|3.5|> and 
conditions (|3.1|) can have only physically meaningful branch points. For simplicity of presentation, let 
us assume that dh < 2a (in the opposite case, a slightly more involved argument still goes through.) 
Using the branch points of g(oj), the only branch points that <f>^ (u>) can have inside (|2.8(l are ±(7r + 
a — Oh). The corresponding branch cuts run along the segments [—n — a — Oh, —n — a + Oh] and 
[rr + a — Oh, ir + a + 0h]- Note that the branch points ±(n + a + 0h) lie outside the physical region (|2.8(l . 
The analogous branch points of $o~(w) can be only ±(7r + a ± icosh - (1/6^)), with the cuts along the 
segments [-7T — a + icosh - (1/6^), — 7r — a — icosh~ (1/0/;)] and [tt + a + icosh~ (l/0h), n + a — 
icosh" 1 {I /0 h )]. 

Indeed, all possible branch points of g{ui) are irn ± Oh (see l|3.7[l). In principle, applying (|3.5fl . 
they could generate many branch points in $f (w) which have no physical interpretation. Let us start 
by showing that —a + Oh is not a branch point: Let us use the fact that $q (cj) is odd to rewrite the 
functional equation i|3.5[l as 

tu(u + a){$([ [a + g(u + a)] - $+[a - ff(w + a)]} + *i 2 (w + a)[$i"(w + 2a) + $]|"(w)] = Qf , (F.l) 
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t 21 (ui + a){^[a + g(u + a)} + $+[a - g{uj + a)}} + t 22 {u + a)[$f (w + 2a) - $f(w)] = Qt ■ (F.2) 
In the vicinity of 0/,, there exist the constants a n such that we have 

oo 

g(u + a) = J2 M" + a ~ Qh) n+l/2 - (F.3) 

n=0 

Since 3>q [a + g>(a; + a)] — $o~[a — .<?(^ + a)] is odd in g, there also exist constants A n such that that 
this function has the expansion 

oo 

4>+[a + g(u + a)} - $+[a - 5 (c^ + «)] = £ A„(w + a - ^)" +1/2 . (F.4) 

n=a 

On the other hand, there exist constants b n such that we can write 

oo 

tn(« + «) = ^ 6 «( w + a - ^)" _1/2 - (F-5) 

n=0 



Thus, the first term in the left hand side of (|F.1|I contains no branch points. The coefficient tm has 
no branch points either. It follows that there are no branch points in Eq. (|F.1|I at all. 

Let us move on to Eq. I|F.2|) . The function $J [a + <?(a; + «)] + <!>([ [a — g(w + a)] is an even function 
of g, and therefore, there exist constants B n such that we can write 



$+[a + g(u> + a)} + $+[a - g{uj + a)} = V B n {u + a- 6 h ) n . (F.6) 



Since the coefficients t 2 i(oj) and t 22 (oj) have no branch points, there are no branch points in Eq. I|F.2(I . 
It follows that the point — a + 9^ which does not have any physical interpretation is not a branch 
point of the function $J"(w). Similarly, it can be shown that the points —a — 9h and a ± 9h are not 
branch points of $f (w). Analogous considerations apply in the antisymmetric case. 

We conclude that the branch points of the Sommerfeld amplitudes of the solution of the original 
problem that lie in the physical region (|2.8|l and therefore give rise to physical waves, lie outside the 
strip Re u el. This means that they do not have to be taken into account in the functional equations 
for 3>f(u)) or by the same token, in the resulting singular integral problem. On the other hand, we 
have no theoretical proof that $>^(uj) that we eventually compute have only physical branch points in 
the physical region i|2.8l) . We can confirm this fact only by carrying out numerical tests. 

Appendix G. The Rayleigh reflection and transmission coefficients. 

When evaluating 1)2.5(1 . a pole 9"^ = a — i/3r of $i(w), with /3r > corresponds to a plane wave 
with the phase factor 

ifcrcos (#— ^iri) „ifcrcos (6 — a) cosh /3r — fcrsin (a — 6) sinh 0b. {(~* 1 \ 

so that its amplitude is exponentially small everywhere except for a small neighborhood of the wedge 
face 9 = a. Thus, we describe a Rayleigh wave incident from infinity along the upper face 6 = a by 
two potentials 

<P^ c (kr, 9) = 47ri0 io e i7fercos(e " e ^ C) , (G.2) 



where </>oo = 1 and <f>iQ — — 2i7RW7j| — 7 2 /(27r — 1); tr = cs/cr with cs being the Rayleigh wave 

speed and 9^ = a — g(i(3s). The reflected wave propagates along the same wedge face as the incident 
but from the tip to infinity and transmitted — along the other face, again away from the tip. They are 
described respectively by 

(t>f{kr,9) = 47rii? rcf ^oe" i7fc, ' cos(6 '" e ^ ) , 

</>f an (fcr, 9) = 4mi? tran io e" i7fcrcos(e+ ^ ) , (G.3) 
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with "scattering angles" 0? R being the complex conjugate of 6^ R C , so that 0q R = a + <?(i/3r) and 
^ir = a + i/^R- Above, i? and i? tran are the Rayleigh reflection and transmission coefficients 

i? rcf = -[R +Tei + R~ Te{ ], i? tran = ^[i? +rcf - iT rof ], (G.4) 

with the symmetric and antisymmetric parts given respectively by 

iJ ±ref = Res[$J; £[(wr) + a], with w R = tt + i/3 R , (G.5) 

so that using the additional angles #o,R = ~ a + 5( w r) an( i $i,R = ~ a + ^r, we have 

As before, the dash denotes the derivative with respect to the argument. 

Appendix H. Adjustable functions and parameters in singular terms. 

As any other numerical code ours relies on a choice of certain options which effect a tradeoff 
between numerical accuracy and either running time or else numerical stability. Apart from the 
relevant grids, these options are 

(i) The adjustable function cr(u>) in f.V.J.5)) In most cases, <r(ui) is chosen to be 

er(w) - -r^r- (H.l) 

sm 75— u> 
2a 

(cf. Budaev and Bogy 1995, Eq. (15)). The choice is convenient, because it simplifies the right 
hand sides of our functional and therefore, integral equations. Also, o~{ui) in i|H.l|) decays at infinity 
reasonably fast. However, for any wedge angle 2a, there exists a critical incident angle 8™ c , such that 
one of the poles of ^ i (w) lies on the boundary of the strip Re en el. For illustration purposes, let it 
be $1 (uj) and let the pole be 80 = tt/2 — a. Then the corresponding term Res ($|f; 8o)cr(ui — 80) has 
one more nonphysical pole, #0 + 2a. In situations like these, another choice of o~(u>) is called for, with 
poles further apart. We have tested 

<r(u;) - -^-. (H.2) 

sm puo 



with various values (3 < n/(2a). However, any a{ui) different from IjH.ll) leads to more cumbersome 
right hand sides of the integral equations and exhibits a slower decay. As a result, the function (|H.2(I 
while increasing the stability of the solution, increases the code run time roughly tenfold. For this 
reason, we abandon i|H.l|) only when 8 lnc is near critical angle. In this region we use l|H.2|) . with 
f3 = ir/(6a). 

(ii) Number of poles. When evaluating the poles of the Sommerfeld amplitudes we do not have to 
restrict ourselves to the strip Re lu e I. The more poles that are utilized in evaluation, the wider the 
domain of analyticity of the corresponding unknowns $j(w) and therefore, the higher the accuracy. On 
the other hand, some nonphysical poles possess residues with large amplitudes and cause numerical 
instability. In the present version of the code, when 8 is away from critical we take into account all 
poles in the strip Re 10 e I, and when 8 is near critical we take into account all poles in the wider strip 
Re u e [n/2 - 2a, tt/2 + 2a] . 

Appendix I. We are grateful to Drs R.K Chapman, D. Gridin and J. Hudson for many useful 
and insightful comments and Professors V. P. Smyshlyaev and A. Gautesen for fruitful discussions 
and suggestions. We would also like to thank Profs. Nazarov and Plamenevskij for invaluable advice 
and relevant references. 
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Appendix J. Nomenclature. 

a(r/) = — i tanh 2??, 



b(v) 



2i sinh ry v 7 2 + sinh 
cosh 2i] 
tanh 2ij 



d{r,) = l . 

XKV) 

e^ (lo) = — tan a sin 2lo + cos 2lo, 
ej" (lo) = tan a cos 2lo + sin 2w, 

, sin w 

ej (w) = tan a cos 2io- 



sin 2w, 



sin lo 



yj^f 2 — COS 2 LO 



Vi 2 - Cl 

ej(w) = tan a sin 2a; — cos 2ui 
g{uo) = cos (7 cos to), 

oo 

nm = Ivp / » 

— OO 

OO 

(W- 1 /)^) = ^tV.P. / C oth^(i -/?)/(*) eft, 
2ai 7 2a 



Wfn= i p /• f(t)x'(t)dt 

nJW 2ai -J S mh^[ X (t)-x(v)}' 



If tanh 2 : 



tanh 2t tanh 2-q I 



^(t-r?) sinh ,£[*(*) -*(„)] J /( ' )dt ' 



<?o"(»7) = ^ 0?) - a(r])Hr+(ri), 

cosh ?7 tanh 2?] 



tfiW) 



tana- 



cosh2?7 2a 

%(v) = ri(v) - b(r])Hr2{v), 



V.P. 



7 cosh rdr 



(1 + 2 7 2 sinh 2 t) sinh ^[r- X (r;)]' 



Qi (v) 



tana- 



cosh?7 



tanh 2r\ 
X' (rj) cosh 2?7 2a\' (jl) 



V.P. 



cosh r dr 



cosh 2r sinh tj^- (t — rj) ' 



rf fa) = - [*o(fl( 2 + ir ?) + a ) ± $ ct(.9( 2 + i? ?) ~ a )\ ~ b ^ [^(^ + ® + iv)± ^t(^ -a + irj) 
rf(v) = ~a(v) [*o(ff(| + ir ?) + a ) T *o(^(| + ir ?) - «)] - [$i"(| + a + ir ?) T *i"(| - a + ir ?) 



rn(w) = -r 22 (w) 



2sin 2w cos lo^J^ 2 — cos 2 w — cos 2 2w 
AM ' 



ri2(w) 
r2i(w) = - 



4cos 2io cos io\J"{ 2 — cos 2 w 

AH ' 

2sin 2o> cos 2u> 



AM 



AM = cos2 2w + 2 sin 2lo cos lo\/~j 2 — cos 2 w, 

X(^) = sinh - (7" sinh(?7)), 

cosh?? 

* W) = T% ■ ,2 =■ 
V7 + smn 77 

Note that the Raylcigh function AM has the purely imaginary root i/3r, with /3r > 0. 
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